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Abstract. We compute the probability of J/ip production from the interaction 
between rj c and the strong magnetic field generated in relativistic heavy ion collisions. 
The computation is first carried out in the heavy quark effective model, in which 
the Ml radiative transition is considered. Then we investigate the transition in the 
framework of non-relativistic heavy hadron chiral perturbation theory and show that 
the polarization of J/ip produced by this process is parallel to the direction of magnetic 
field and thus perpendicular to the reaction plane. The transition probability obtained 
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1. Introduction 

The production of heavy quarkonia in relativistic heavy ion collision has been widely 
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in . 

CN . studied in consideration of their possible role as probes of the quark-gluon plasma (QGP) 
and the QCD deconfinement transition. J/ip is the most promising candidate among 
qq • the charmonium probes because of its relatively large production cross section and its 
large branching ratios for leptonic decays. Various competing mechanisms contribute to 
the production and the suppression of J/ip in nucleus- nucleus collisions. The so-called 
"anomalous" suppression of J/ip was first observed in experiments at the CERN-Super 
Proton Synchrotron (SPS) by the NA38 collaboration [1], then confirmed by the NA50 
and NA60 collaborations [2,3], and later observed in experiments at the Relativistic 
Heavy Ion Collider (RHIC) by the PHENIX collaboration [4]. 

Theoretically, charmonium states are successfully described as nonrelativistic bound 
states of a cc pair with a binding potential incorporating one-gluon exchange and a 
linear confining potential. Matsui and Satz [5] proposed that the color screening of 
the quark- ant iquark potential in a quark-gluon plasma will hinder the formation of 
charmonia. Their prediction of the suppression of charmonium production preceded the 
experimental observations; however, the measurements did not confirm the anticipated 
suppression pattern. In order to understand the experimental data in detail, various 
other mechanisms for charmonium suppression have been investigated, so far without 
a generally accepted, comprehensive solution. In parallel, the color screening effect has 
been studied in several different approaches. Recent results on color screening based on 
calculations in lattice QCD can be found in [6]. 

Besides the screening effect, the suppression may be caused by the kinetic 
dissociation of cc pairs. Due to finite formation time of charmonia, while a cc pair is 
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produced in the collision, the pre-charmonium state may interact with the hard gluons 
when traversing the medium and dissociate [7,8]. In the hadronic phase, the cc pair 
may combine with light quarks to form a pair of heavy mesons [9]. Moreover, the 
production of J / ip can potentially be influenced in the late-phase hadronic environment 
by scattering with other hadrons and dissociation into a DD pair [10]. These effects 
will also contribute to the suppression. 

In relativistic heavy ion collisions, cc pairs are created in the hard scattering of two 
gluons, and different charmonium states may be formed when the resonant energy is 
reached and the correct quantum numbers are established. For example, according to 
Yang's theorem [11], t] c and J/ip can be produced through gg — > rj c and gg —> J/ijjg 
respectively. Nevertheless, not all charmonia are directly produced in the collision. As 
the energy of collision increases, the number of cc pairs created by hard scattering grows. 
The recombination of these charm and anti-charm quarks from different binary nucleon 
collisions may cause a considerable enhancement of the charmonium production. Based 
on the hydrodynamic simulation of the central collision under RHIC condition, between 
30% and 50% of the J/i/j yield may come from the recombinant production [12,13]. 
For instance, through the measurement of the e + e~ decay mode, the decay chain of 
Xciip') ~~ >* J 1^1 ~~ >* e + e~7 can be reconstructed in p+p collision [14]. Furthermore, the 
suppression of excited charmonia may also result in the reduced production of J / tp and 
other low energy states in the charmonium spectrum due to their greater fragility in the 
medium. 

Motivated by the J/ip production through charmonium decays, we here investigate 
the J/ijj production through the radiative Ml transition of r] c induced by the strong 
magnetic field generated in the ultra-relativistic A+A collision. The ability of such 
transient fields to mix pseudoscalar and vector mesons in the light quark sector was 
recently discussed by Asakawa et al. [15] in connection with the chiral magnetic effect 
[16]. The same mechanism also induces a magnetic dipole transition between the 
pseudoscalar and vector states of heavy quarkonia. Here, we first compute the transition 
probability t] c + B — > J/tp, where B denotes the magnetic field, by considering the 
magnetic dipole interaction in the quark model. Then we confirm the computation by 
utilizing non-relativistic heavy hadron chiral perturbation theory (NRHH\TT [17]) to 
the leading order. Finally, we estimate the transition probability in a Au+Au collision 
at RHIC. 

2. Radiative Ml transition 

To analyze the transition of J/ip from the interaction between r] c and magnetic field 
generated in the collision, we consider the corona scenario which applies to the hadronic 
matter remaining near the surface of the colliding nuclei [15]. An r] c meson produced 
near the surface can still interact with the magnetic field generated during the early 
time when the magnetic field is very strong. We choose the x — y plane as the reaction 
plane and assume that the beam is oriented along the x direction. We also assume the 
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nuclear charge Ze to be distributed uniformly within the nuclear radius R. In peripheral 
collisions, when the impact parameter b « R, the magnetic field at x = on the surface 
can be approximated as 

Zejvb A . . 

where 7 = 1/ a/1 — v 2 is the Lorentz contraction factor. The details of approximation 
can be found in the Appendix of [15]. Since the collision is ultrarelativistic, the two 
ions can be visualized as two colliding pancakes due to the Lorentz contraction in the x 
direction. Therefore, we expect the region to be in the vicinity of x = 0. 

We now outline the computation of transition probability in the quark model. Since 
rj c and J/ip are pseudoscalar and vector mesons, respectively, conservation of parity 
and charge conjugation symmetry dictates that the transition will be governed by the 
magnetic dipole interaction. By using the time dependent perturbation theory, the 
transition probability amplitude can be written as, 

a{t) = -if e^ 1 ' {f\HUt')\i)dt\ (2) 

where E ^ = E f — Ei is the energy difference of the initial and final states. In our case, 
(f\ = (J/i>l \i) = \Vc) and H int (t) = -fj, c -B(t), where fx c = ^(S c -S g ) is the magnetic 
dipole moment of the cc pair. Here e c = 2/3 and m c denote the charge and mass of the 
c-quark, respectively, and S c ( g ) denotes the spin of the c-(c-) quark. We neglect strong 
interaction corrections to the g-factor of c-quark. By taking the limit t — > 00, we obtain 
the total transition amplitude in the quark model approach as 



fl 4vr 3m c UI cz czl ' loo (R 2 + v^H' 2 )^ 



Zejvb 2e 2E fi K 1 (^ 



(3) 



An 3m c v 2/ y 2 R 
where K\ is the modified Bessel function of second kind and 

(f\S cz - S- cz \i) = (J f ,M f \S cz - S 5z \J t ,M z ) = 1 (4) 

for Jj = 0, Jf — 1 and Mj = M/ = 0. The transition probability is given by 

QM12 = 16vra Z 2 ai 2 v 2 b 2 ( 2E fi K 1 2 

1 9m 2 c 4tt 1 v 2 ~f 2 R ) ' U 

The prefactor 167ra/(9m 2 ) represents the strength of the coupling between the 
charmonium spin and the magnetic field. A similar result will be expected in the 
effective field theory (EFT) calculation, which we perform next, while the prefactor 
will be replaced by a coupling constant which encodes the non-perturbative physics of 
the charmonium bound state. 



J/ip Production by Magnetic Excitation of r) c 



4 



3. Calculation in the effective field theory 

To calculate the transition probability in the formalism of EFT, we introduce an effective 
Lagrangian that describes the interaction between rj c , J/ip and the magnetic field. 
Since the hyperfine splitting of r] c and J/ip is about 100 MeV, the transition can 
be described non-relativistically in the rest frame of the J/ip. We thus make use of 
the framework of non-relativistic heavy hadron chiral perturbation theory (NRHHxPT 
[17]). By preserving parity, charge conjugation, and chiral symmetry, the leading order 
interaction term can be written as, 

C mt = ^Tr[JB-aJ% (6) 

where p is a coupling constant and J = ip ■ a + rj c is the rj c — J ftp chiral multiplet. Due to 
the presence of the Pauli matrix, this term can break the heavy quark spin symmetry. 
Then the transition amplitude can be written as 

T fi = (J/i!,B cx \r] c ,B cx ) 

= ig' J alpfdpi 0}/^(p/) <P Vc (~Pi) 

d 4 x (J/^, J B cx |^(x)?7 c (x) J B At (x)|?7 c , J B cx } 

= ig' J dpfdpi 0}/^(p/)<MPi) 
x (7/^,^1^(0)^(0)1^,^) 
x J d 4 xe^f- p ^B cx (x)^ 

= ig' J dp f dpi 0} /0 (p/) (j)r, c (pi)B ex (p f - Pi) ■ e JA/ \ (7) 

where dpf = d 3 p j / [(2tt) 3 ^J2E f] is the relativistic phase space measure, and 4>j/^(Pf) 
and (pnXVi) represent wave packets of the J/ip and r\ c . By taking the plane- wave 
approximation, <pj/^(pf) and (j) Vc (Pi) will be simply identity. The subscript /i denotes 
the direction of polarization, e J '^ is the polarization vector of J/ip, and B eyi is the four- 
dimensional Fourier transform of the external magnetic field. Finally, g' = 2 y jmj/^m, nc p 
takes into account the normalization of the non-relativistic heavy fields. The unknown 
parameter p can be extracted from analyzing the decay width of J/ip — > 'yrjc in quark 
model [18]. 

In order account for the Lorentz contraction of the magnetic field (1) in the re- 
direction, we now keep its full x-dependence. For simplicity, we consider the production 
at y = near the nuclear surface. The z-component of the magnetic field contributed 
by a single nucleus can be written as 

B&(t x) = ^ (8) 

2 l ' ] 87t( j R 2 + 7 2 (x T ^) 2 ) 3 / 2 ' K ) 

where the minus (plus) sign is for a right- (left-)moving nucleus. We firstly evaluate the 
contribution for the right-moving nucleus. By taking v = 1 and choosing the light-cone 
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(E fi +p fi )R 



x S\ Pf± - PiX )(E fi+Pfi ) ^±-^ L ( 9 ) 

where p^ = pf x — pi x and 8 (Efi — pfi) can be rewritten as 

(Ej-Ptef +pj ± + m 2 f / p 2 f± + mj - {Ej-p lx ) 2 \ 
2(E i -p ix y 6 { Pfx 2(E i -p ix ) )■ 

Since mj/^p m m Vc 3> Pf,Pi non-relativistically, we will make the following approximation 
to simplify the expression: 

/ p 2 f± + m) - (Ei - p ix ) 2 \ ( m 2 f - m 2 + 2p ix mi \ 
°\Pfx TzfB \ «<MP/«- 



2{E i -p ix ) J V 2m, ; 

~ S(Pfx~Pix- E fi) (11) 

and 

{Ei - p ix f + p 2 fl _ +m 2 f ^ m 2 j/ip + m 2 Vc . (12) 
Inserting B cx (pf — p.j) into (7), we obtain the transition probability amplitude, 
(+) ,Zejb f d Pi fm 2 j^ + m 2 \ 



a f ~ W / ~r^= <Pj/f{Pf) VvcKPi 



to J ^2E f TJIWxrJ ' T " cx ^' V 2m l J 

x2E fi T^i \pf=iPix+E fi , Pi± )- (13) 

Repeating the calculation for the minus sign in (8), we find a\ + ^ — a\ . If we only 
keep the leading-order term in the expansion of Efi, the full nonrelativistic transition 
probability can be expressed as, 



\a fi 



NR|2 _ 9' 2 (Ze lb y ( 2E ft K 1 (^ 
AE t E f \ 4tt J \ -f 2 R 

Z 2 ai 2 b 2 hEfjKj^ R) 
j/4> 4tt ^ ^R 



which is equal to the transition probability one finds by just using the magnetic field 
expression (1). This correspondence reflects the fact that produced J ftp mesons are 
highly localized at x = 0. 

In the previous calculation, we analyzed the transition in the rest frame of J/ip. 
Nevertheless, the r\ c directly produced in the collision may move relativistically. Under 
this circumstance, the NRHHxPT has to be promoted to its relativistic version[19] and 
the charmonium field must be rewritten as 

J=^(1+^)[^7 M -^75]^(1-^), (15) 

where is the four-velocity of the charmonium multiplet J. However, even in the 
boosted frame, the structure of the interaction term remains unchanged. Although the 
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transition matrix element only depends on the energy of photon, the transition rate can 
vary due to the kinematics. When r\ c carries large transverse momentum, and we go 
into the longitudinally co-moving frame of the charmonium states, we have the relation 

Ip/j.1 ~ IPiJ > 'mj/ il ,m Vc ^>p fx ,p ix (16) 

along with the delta function which imposes the conservation of transverse momentum 
in (9). Thus, we can approximately write 

E fl ^(m 2 j / ^-m 2 ric )/(2p i± ), 

{Ei - Plx ) 2 « p 2 ± + m 2 - 2p ixPi± . (17) 
Using these approximations, the relativistic transition probability becomes: 

g' 2 fZejbV (2E fi K 1 (^] 



\ a fi I 




Compared with the nonrelativistic result (14), the relativistic result (18) contains a 
factor p 2 ± in the denominator, which causes a significant suppression when the transverse 
momentum is much greater than the charmonium mass. 



4. Numerical results and conclusion 



Now, let us find some numerical results for Au+Au collisions at the highest RHIC energy. 
By taking Z = 79, b = R = 7 fm, m c = 1.4 GeV, v — 1 and 7 = 100, and using (5), we 
obtain the transition probability |ctjj M | 2 = 2.37 x 10~ 4 in the quark model. To derive 
the numerical result for the EFT calculation, we need to determine the numerical value 
of p. By computing the decay width of J/ip — > 7% in the framework of NRHHxPT and 
comparing the result with that derived in the quark model [18], we find p = 0.147 GeV" 1 
and thus |ajj R | 2 = 2.37 x 10 -4 . If the relativistic correction is taken into account, we find 
p = 0.137 GeV -1 and |ajf| 2 = 2.06 x 10~ 4 . For the transition with large energy of the 
photon, the probability amplitude can also be computed in EFT framework by utilizing 
potential non-relativistic QCD(pNRQCD)[20], while the result at leading order would be 
tantamount to that in the quark model. We thus conclude that the leading-order result 
for the magnetically induced transition probability from n c to J/ip is about 2 x 10 -4 . 
Unfortunately, this contributi on to J/^-production is hardly observable experimentally. 

We now consider some possible sources of corrections to this result. At the next-to- 
leading order, the coupling constant p can be affected by the contribution of intermediate 
hadronic loops, which is indicated in the recent study of charmonium decays [21,22]. 
However, when matching to the experimental data of J/ip — > / ~fr] c decay, the correction 
will lead to further suppression on p, which makes the transition probability even smaller. 
In pNRQCD, the correction from the higher derivatives will result in the suppression of 
transition matrix element as well [20]. According to (14), the transition probability will 
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be independent of 7 in the ultra-relativistic limit (7 — > 00) since K\{x) ~ x _1 + 0(x) as 
x — > 0. Increasing the collision energy further (e. g. to LHC energies) will thus not change 
the result much. On the other hand, the result (18) for large transverse momentum will 
give a more substantial suppression. We thus conclude that our quantitative estimate 
represents a robust upper bound. 

In summary, we analyzed the contribution to J/ip production in relativistic heavy 
ion collisions through the magnetic dipole transition of t] c due to the interaction with the 
strong magnetic field generated in the collision. We evaluated the transition probability 
in the framework of the quark model and in the EFT approach. The EFT approach 
shows that the polarization of the produced J/ip will be parallel to the direction of 
the magnetic field. Both approaches lead to comparable results, in which the transition 
probability is about 2 x 10 -4 , which is probably too small to be observed experimentally. 
We also considered the case when the r] c and J/ip carry a large transverse momentum 
and found that the result is further suppressed. 

This work was supported in part by a grant from the U.S. Department of Energy 
(DE-FG02-05ER41367). 
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